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Abstract. A Markov evolution of a system of point particles in R d is de- 
scribed at micro-and mesoscopic levels. The particles reproduce themselves at 
distant points (dispersal) and die, independently and under the influence of 
each other (competition). The microscopic description is based on an infinite 
chain of equations for correlation functions, similar to the BBGKY hierarchy 
used in the Hamiltonian dynamics of continuum particle systems. The meso- 
scopic description is based on a Vlasov-type kinetic equation for the particle's 
density obtained from the mentioned chain via a scaling procedure. The main 
conclusion of the microscopic theory is that the competition can prevent the 
system from clustering, which makes its description in terms of densities rea- 
sonable. A possible homogenization of the solutions to the kinetic equation in 
the long-time limit is also discussed. 



1. Introduction 

Phenomenological physical theories of substances like gases, liquids, etc, are 
based on macroscopic observations and thus consider these substences as continu- 
ous media. However, already in the time of L. Boltzmann and M. Smoluchowski 
it was understood the importance of the microscopic structure of such and other 
similar substances. Since then a challenging problem of theoretical and mathemat- 
ical physics has been to derive the rules of the collective behavior of large systems 
of interacting particles from their microscopic theory, based on the so called 'first 
principles'. Achievements in this direction laid the fundamentals of modern statisti- 
cal physics; see, e.g., [10]. In the Hamiltonian mechanics, the motion of N physical 
particles in R d is described by a system of 2dN differential equations, subject to 
initial conditions. For N 3> 1, the abundance of equations, and hence of the initial 
conditions, makes the point-wise description meaningless since no observation could 
indicate at which point of the phase space the system actually is. Moreover, the 
point-wise description would be 'too detailed' for understanding the collective be- 
havior of the system. An alternative can be the statistical approach which provides 
the possibility to link micro- and macroscopic descriptions to each other. In this 
approach, one deals with probabilities with which such points lie in given subsets 
of the phase space. Then instead of 'point- wise evolution' one studies the evolution 
of probability measures, considered now as the states of the system. In 1| , N. N. 
Bogoliubov suggested another approach, which later became popular, especially in 
the physical literature. Here, the evolution of probability measures is described 
indirectly as the evolution of the so called moment or correlation functions. The 
latter evolution is obtained from a hierarchy of equations, called now Bogoliubov 
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or BBGKY hierarchy or chain, that couples correlation functions of different order 
with each other, see |10j . The description at this level is microscopic since one deals 
with coordinates of individual particles; cf. the Introduction in [TB]. More coarse- 
grained levels are meso- and macroscopic ones. They are attained by choosing 
appropriate scales for space and time, see [THUS]- Of course, certain details of the 
system's behavior are lost due to such coarse-graining. Kinetic equations provide a 
space-dependent mean-field-like approximate^ description of the evolution of large 
particle systems. Such equations are deduced from the BBGKY chain by means of 
various types of scaling, cf. Section 6 in [TU] and also jTSJ . 

Along with physics, in modern life sciences one also deals with systems con- 
sisting of large number of interacting entities distributed in a continuous habitat 
and evolving in time. Their collective behavior is observed at a macro-scale, and 
thus the corresponding mathematical theories which explain this behavior tradi- 
tionally describe their dynamics by means of phenomenologically (or heuristically) 
deduced nonlinear equations, mostly differential or integro-differential, involving 
macroscopic characteristics like density, mobility, etc. However, this kind of macro- 
scopic phenomenology may often be insufficient also in life sciences, as they can col- 
lect massive experimental data of high precision characterizing individual behavior 
of constituting entities. For instance, technical advances in animal tracking and 
geno typing allow for obtaining rich individual-level data on population dynamics, 
population genetics, and evolutionary biology. However, mathematical methods 
are absent that would enable one to derive the population-level consequences of 
such individual-level observations. From this and many other examples one may 
see how important is to elaborate a general mathematical theory in which the rules 
of the collective behavior of a large number of interacting entities can be deduced 
from the information about their individual behavior and the interactions. This 
article is aimed at announcing some of our recent achievements in this field based 
on the concepts of statistical physics mentioned above. The main principles which 
we follow in our study are: 

• The objects are large systems of entities distributed over a continuous habi- 
tat (part of M. d , d > 1) and evolving in continuous time. 

• The phase space of the model is the collection r of possibly infinite sets of 
identical points - configurations 7 C M. d . 

• In the microscopic theory, the states are probability measures on the phase 
space r, the evolution of which is assumed to be Markovian, described by 
Kolmogorov-type equations. The main acts of the evolution are particle's 
birth, death, and spatial motion - jumps or diffusion. 

• The evolution of states is described indirectly as the evolution of the corre- 
sponding correlation functions obeying equations similar to BBGKY chains. 

• The mesoscopic theory is based on kinetic equations, deduced from the 
equations for correlation functions by means of scaling procedures, in which 
one scales interactions and densities, whereas time is left unsealed. 

2. The setup 

The object of our study in this paper is a large population of motionless entities 
(e.g., perennial plants) which reproduce themselves, compete and die at random. 



In the physical literature, it is rather called random phase approximation. 
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The reproduction (dispersal) consists in random independent sending by an entity 
an offspring to a distant point, which immediately after that becomes a member of 
the population. Each entity can die independently with a constant rate, as well as 
under the influence of other population members. The latter event is interpreted 
as competition-caused. In theoretical biology such models were introduced and 
studied, e.g., in (2j[3j[l6]. A more detailed bibliography can be found in [6]. 

As was suggested already in [2], the right mathematical context for studying 
models of this kind is the theory of random point fields in R d ; cf. also page 1311 
in [15] . Herein, populations appear as particle configurations constituting the set 

(1) r := {7 C R d : I7 n K\ < 00 for any bounded K C R d }, 

where, for a set A, by |^4| we denote the number of points in A. For our model, T 
is a phase space. Its metric properties, as well as the main aspects of the analysis 
on such spaces, were studied in [11]. Along with finite ones the set T contains also 
infinite configurations, which allows for describing 'bulk' properties of a large finite 
system ignoring boundary and size effectfl If the initial configuration 70 is fixed, 
the evolution might be described as a map t > 74 £ T, which in view of the random 
character of the events mentioned above ought to be a random process. However, at 
least so far, for the model considered here this way can be realized only if 70 is finite. 
Hence, to describe mentioned bulk properties one ought to go beyond the point- 
wise setting, as it has been done in the statistical physics of interacting particle 
systems. As mentioned above, in the statistical description states of the system 
are probability measures on T. To characterize them one employs observables - 
appropriate functions F : T — > BL As in the Heisenberg approach in quantum 
physics, the evolution of states can be described via the evolution of observables. In 
view of the assumed Markov property, this evolution is obtained from the equation 

(2) -F t = LF t , F t \ t=0 = F , t>0, 
at 

in which the operator L determines the model. In our case, it is 

(3) (LF)( 7 ) = £[m + £r(s >7 \ a: )][*'(7\ a; )- J F(7)] 

+ f E + (y, 1 )[F( 1 Uy)-F( 7 )]dy, 

JR d 

where 

(4) E ± (x, 1 ):=J2"±(x-y)- 

ye-y 

Here a± > are suitable functions, see (fTU|) below. As given in ((3]), L is a typical 
'birth-and-death' generator in which the first term corresponds to the death of the 
particle located at x occurring (a) independently with rate m > 0, and (b) under 
the influence of the other particles in 7 with rate E~~(x,j \x) > 0. Here and in 
the sequel in the corresponding context, we treat each x € M. d also as a single-point 
configuration {x}. That is, if x belongs to 7 (resp. y does not), by 7 \ x (resp. 
7 U y) we mean the configuration which is obtained from 7 by removing x (resp. by 
adding y). The second term in ([3]) describes the birth of a particle at y g M. d given 



A discussion on how infinite systems provide approximations for large finite systems can be 
found in [4] [5]. 
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by the whole configuration 7 with rate E + (y, r y) > 0. The model described by ([3]) 
is called a spatial logistic model (SLM). A particular case of SLM is the continuous 
contact model [TU [T3] where a_ = 0, and hence the competition is absent. 

The sums in ([3]) and ([4]) may contain infinite number of summands and hence may 
diverge. This means that the direct study of the corresponding equation (J2j) for our 
model is rather problematic. Following the way suggested by N. N. Bogoliubov, the 
evolution of states will be studied via the evolution of the corresponding correlation 
functions ko > k t . By a correlation function k we mean an infinite collection of 
symmetric positive functions k^ : (R d )™ — > R, n £ No, which determine the state 
in a certain way, see [Sj [Tl [Sj [9] for more detail. Such k can also be viewed as infinite- 
dimensional (Fock-type) vectors. Note that k^ = 1 and k^ is the particle density. 
For the Poisson measure with intensity h > 0, one has k^ n \x\, . . . , x n ) — >c n , 
n € No- We say that a correlation function k, and hence the corresponding state, 
are sub-Poissonian if 

(5) k^\xi,...,X n )<C n , 

which holds for some C > 0, all n £ No, and (Lebesgue) almost all (xx, . . . , x n ). 
Then a sub-Poissonian state is, in a way, similar to the Poissonian state in which 
the particles are independently scattered over R d . At the same time, the increase 
of fcW with n as nl (see (JSJ below), may correspond to the formation of clusters. 
The evolution fco i— > k t is obtained from the equation 

(6) ^_k t =L A k t , k t \ t=0 = k Q , 

which, in fact, is a chain of equations for particular k\ , analogous to the Bogoli- 
ubov hierarchy for physical particles. The generator L A is constructed from L as 
in ([3]) by a certain procedure, see |9]. In particular, the first equations in the 
chain © have the following form: dk^ /dt = and 

,7, 4<»(*) --**«(.)-/ M-»)*f (*»)*+/ 

JR d JR d 

The RHS of the equation with dk^'/dt contains k[ n ^ with n — 1,2,3, etc. Theo- 
retical biologists try to solve chains like (j6)) by decoupling them; cf. [14] . In the 
simplest version, one sets 

k^(x,y) ~k w (x)kM{y), 

which is equivalent to neglecting spatial pair correlations. This turns ([7]) into the 
following nonlinear equation 

which, in fact, is a kinetic equation for our model. Its eventual solution provides 
a mean-field-like approximation of the evolution of the density. Note that the 
question of whether the evolution of states preserves sub-Poissonicity, and hence 
clusters do not occur, can be answered only by studying the whole chain in (p]). 
For the contact model, which one obtains from ([3]) by setting a_ = 0, it is known 
[7] that 

(9) const • nl c" < fc^ (xi, . ..,x„)< const • nl C™, 



STOCHASTIC EVOLUTION OF A CONTINUUM PARTICLE SYSTEM 



5 



where the left-hand inequality holds if all x% belong to a ball of small enough radius. 
Hence, in spite of the fact that Ct — >• as t — » +00 if the mortality dominates the 
dispersal, kt are definitely not sub-Poissonian if a_ =0. 
In this article we address the following questions: 

(1) Can the competition in §S§ provide that the problem © has a solution that 
obeys ©? 

(2) Can one deduce © rigorously; e.g., by a kind of scaling procedure, which 
would allow for error control? 

(3) What can be said of solutions to ((SJ in that case? 

The answers are given in the form of theorems stated below. The proof of Theorem 
[T] can be found in [6] ; the remaining theorems will be proven in our work under 
preparation. 

3. The results 

We assume that the competition and dispersal rates are additive, cf. Q, and 
the corresponding kernels a± are symmetric, essentially bounded, and integrable. 
That is, 

(10) a± G L 1 (M d ) fl L°°(M. d ), a±(x) = a±(-x) > 0, 
and thus we set 

(11) (a±) = / a±(x)dx, \\ a ±\\ = ess sup a± (x) . 

Jw d xes. d 

Let Tq be the subset of T as in ([I]) consisting of finite configurations only. For 
7 G To, by I7I we denote the number of points in 7. One observes that, for 7 G To, 
the sums in ((J]) and in the expressions below are well-defined. Set 

s ± ( 7 ) = ^i? ± (x,7\x) = ^ a ±( x ~y)> 7er . 

2:67 yei\x 

Then the operator in ^ acts as follows 

(12) (L A fc)( 7 ) = -(m| 7 | + E-( 1 ))k( 1 ) + Y,E + (x, 1 \x)k( 1 \x) 

- / B-tomW* I £« + <*-„«7\zU„)*, 

where 7 G To. For a one-point 7 = {x}, one readily derives from (|12p the right-hand 
side of 0. 

So far, f|12[) is just a formal expression. Then our next aim is to introduce a 
Banach space, in which L A can be defined as a linear operator, and hence (j6|) can 
be studied. For k : Tq — >• R, let us assume that each k^ : —5- R is essentially 
bounded, and then, for a G R, set 

(13) ||fc||„ = sup e an g n (fe), 

neN 

where go(fe) = |^ I, and 

q n (k) := esssup |fc (Tl) (xi, . . . , x n )\, n G N. 

(a;i,...,a;„)e(R £i )™ 
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Afterwards, we set 

(14) K a = {k :r Q ->M : \\k\\ a < oo}. 

That is, K. a is an infinite dimensional real vector space equipped with norm (|13p. 
complete with respect to this norm. For a" < a', we have ||fc|| Q » < ||fc|U'; and 
hence, 

(15) K a > C lC a », for a" < a'. 

That is, the bigger is a, the 'smaller' is the corresponding space (TT4"1) . For k g JC a , by 
(|13p we readily have q n {k) < ||fe[|ce exp(— an), which yields that /c is sub-Poissonian; 
cf. JS|). Thus, for each a£l, the space (fl"4")l contains only sub-Poissonian fc, and 
each such a A; is contained in some K, a . Therefore, the collection of spaces (jT4"]) with 
all possible a£l, cf. (1151) . provides an appropriate framework for describing the 
evolution of sub-Poissonian correlation functions. 

Now let us turn to defining the operator (fT12j) . First of all, we note that, for 
k 6 JC a , L A k need not be in the same JC a , which means that L A , as a map from 
JC a into itself, cannot be defined on the whole /C Q . On the other hand, it is possible 
to obtain the following estimate, see Section 4 in [5J, 

||£ A fc||a < 4>{ a ' — a )\\k\\a', a' > a, 

where ip : (0, +oo) — > (0, +oo) is a certain function such that ip(r) — > +oo as 
t — > + . That is L A acts as a bounded operator to JC a from any smaller K. a ' . Then 
the domain of L A in K, a is defined to be the set of all those elements of this space 
which belong to smaller spaces IC a > . That is 

a' :q' >q 

Now we can make precise what do we mean by solving ©. Namely, for < T < 
+oo, a map [0, T) 3 t i— > fcj G /C Q is said to be a classical solutior^ to the problem 
([B]) in space K, a on time interval [0, T) if: (a) it is continuously diffcrentiablc on 
(0,T); (b) k t S P Q for all £ £ (0,T); (c) it obeys both the equation and the initial 
condition in (JSJ. 

Let us take some a* 6 R, and then any a* > a*, and then set 

(16) r* 



<a+) + (a_)e-«« ' 

where (a±) are the same as in (jlip . Our first result is given in the next statement. 
Theorem 1. Suppose that there exists 8 > such that, for almost all x € M. d , 

(17) a + {x) < 6a-{x). 
Then, for each ci'cR such that 

(18) 8e a * < 1, 

and any a* < a* , the problem with ko € /C Q * /las a unique classical solution in 
/C a „ on iime interval [0, T*), w?z</i T* given in H16\) . 



'In infinite dimensional Banach spaces, there can be solutions in weaker senses. 
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Let us make some comments on this result. First of all we note that the solution 
to (|6|) is sub-Poissonian, as hence clusters do not occur, if the initial fco shares this 
property. Second, the solution lies in a 'bigger' space than fco does. The choice of 
this initial space is restricted by (I18|) . For a fixed a* , T* has a unique maximum as 
a function of a*, that means that there exists an optimal choice of a*, for which 
the time interval [0,T* ) has maximal length. In (|17j) we, in fact, assume that the 
dispersal is dominated by the competition. The essence of this condition can be 
formulated as follows: each particle can 'kill' its offspring and can be 'killed' by it. 

Now let us turn to the mesoscopic theory, which we base on a Vlasov-type scaling; 
cf. Section 6 in [5] . According to the general scheme developed in [8] , this amounts 
to passing to the scale at which the microscopic density gets of order £ _1 and the 
interaction of order e, where e € (0, 1] is a scaling parameter. This means that the 
initial correlation function fco (7) gets of order e - ' 7 '; and hence, one assumes that 

(19) £ l7 'fco(7) ->r (7), as e -> 0+ 7^0, 

where r is a correlation function. Additionally, one assumes that also fc t for t > 0, 
have the same divergence, and hence should be 'renormalized' as in (TT^) . The 
interaction in L A is represented by the competition kernel a_ which thus should be 
replaced by ea_ (weak interaction limit). Let L A be the operator which we obtain 
by such a replacement. Then we set, cf. (|19|). 



(20) fcg en ( 7 ) = ^lfc t ( 7 ), (^ A rc „fc)(7) = e hl (L? k)(j). 

Thereafter, the problem ([6]) turns into the following 

(9-n —h {e) - T A h {e) h {£) I - fc (e) 

dt * Ten ~ e > ren *> rcn ' fi, t,ronl*=0 — ^Ojren" 

The expression for L^ lcn can be calculated from ([2"0]) explicitly, which yields it in 
the form 

(22) L* Ien = V + eB, 

where, cf. (|T2|) . 



(23) (Vk)(j) = -mk(y)- E~(x, -y)k(-y U x)da 

+ / y^a+ix - y)k(j\xUy)dy, 



and 



££7 



(Sfc)( 7 ) = -E- ( 7 )fc(7) + J2 E+ ( X ^\ x )Kl \ x). 



Note that L^ rcn | e= i = L A . For each e € (0, 1], (|2"Tj) has a unique classical solution 
in the same space and on the same time interval as ([5]) does. This can be proven 
exactly as Theorem [TJ For r$ as in (119)) and V as in ([2"3"]l . let us consider 

(24) ^-r t = Vr t , r t | t=0 =r . 

dt 

For this problem, one can probe an analog of Theorem [TJ even without assuming 
(IT71) because the latter condition relates to the operator B, cf. 
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Theorem 2. Let a*, a*, and be as in Theorem^ and let k[ e } cn (resp. r t ) 
be solutions to [2b]) (resp. \2J$ ) in JC a _ t on time interval [0,T*) ; with the initial 
conditions k[ e ^ cn \ t= o — r t \ t=Q = ?*o G IC a » . Then, for every positive T < T*, it 
follows that 

sup \\k\ e } en - r t \\ a , -> 0, as e -> + . 
te[o,T] 

The equation in is, in fact, a chain of equations for r^ 71 ' , n 6 No- It is called a 
Vlasov hierarchy. Theorem states that, in the scaling limit e — > + , the solutions 
to converge in /C Q , to the solution to the Vlasov hierarchy, uniformy on closed 
subintervals of [0,T*). The reason to pass from fc t = k[ r en |e=i to r t = k[ r en |e=o can 
be seen from the following arguments. Suppose that r in (|24[) is such that rg ' = 1 
and 

r^\xi, ...,x n ) = go (a;i) • • • £o(#n) 5 n € N, 

for some positive go £ In other words, we assume that ro is the correlation 

function of a heterogeneous Poisson measure, with density function q$. Let us now 
seek the solution to (f!M)) also as the product of certain unknown g t . The peculiarity 
of V as in ([25]) is that it is possible to do, and that g t can be found from the 
following equation 

(25) ^-Qt(x) = -mgt(x) - Qt(x) f a-(x - y)g t (y)dy + [ a + (x - y)g t (y)dy, 
at j R d j K d 

with the initial condition Qt\t=o — Qo, which is exactly the kinetic equation ([Sjl. In 
particular, the passage from e = 1 to e = has led us from (J7J to ©. Since we 
know that the Vlasov hierarchy (|24l) has a unique solution, then this solution has 
the form r[ n \xi, . . . , x n ) — gt(xi) ■ ■ ■ gt(x n ), n G N, whenever (|25|) has a unique 
solution in L°°(M. d ) on time interval [0, T»), which obeys the condition 



||0t||L~(K<i) < e~ 

The latter is to guarantee that r t lies in /C Q „; cf. (fT3")) and (HU). 



Our next problem is to study the solvability of (|25|) . For some technical reasons, 
it is more convenient to solve it in the space Cb(R d ) of all bounded continuous 
functions (j> : M. d — > R. We equip this space with the norm 

||0|| = sup \(f>{x)\, 

x£R d 

which turns it into a Banach space, that can be isomorphically embedded into 
L°°(R d ). We say that tf> G C h (R d ) is positive if <j>{x) > for all x € K d . Our first 
result in this domain is the following statement. 

Theorem 3. For each positive qq G Cb(M d ), the problem \25\) with the initial 
condition Qt\t=o = Qo has a unique classical positive solution g t G Cb(K d ) on time 
interval [0, +oo). 

Now let go G Cb(R d ) be constant as a function of x. Then the solution as 
in Theorem [3] will also be independent of x and can be obtained explicitly. Set 
gt{x) = u t . Then u t has to solve the following 

(26) ^j-u t = ((a+) - m)ut - (a-)u^, 
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see (jllj) . which is a Bernoulli equation. For m > (a+), its solution decays to zero 
exponentially as t — > +00. For m = (a+), the solution is ut = uq/(1 + (a-)uot), 
and hence also decays to zero as t — > +00. For m < (ct+), we set 

(a+) — m 
q= <«-> ' 
In this case the solution to (|2l)|) has the form 

u q 



u Q + (q - u j cxp(-g(a_)<) 
which, in particular, means that u t q as t — > +00. 

Theorem 4. Suppose that q > and, /or almost all x G R d , £/ie following holds 

"+{x) ( m \ a-(x) 



(28) -fY > 1 

(a+) V («+)/ («-) 

Lei a/so the initial condition go G Cb(M d ) to the problem (25)) obeys 

< S < q (x) < q, 

which holds for some S and all x € M. d . Then, for each x £ M. d and t > 0, the 
solution as in Theorem^ obeys the bounds u t < Qt(%) < q, where Ut is given in 
\2T§ with uq = 5. Hence Qt{x) — > q point-wise as t —> +00. 

The result just obtained can be interpreted as an asymptotic homogenization 
of the density. Let us add some comments to (f28|) . If a+(x) = 9d-(x), for some 
9 > and almost all x, cf. ([17]). then flU} holds for all m G [0, (a+)). Note that 
the case of a+(x) = a~(x) was studied in [T7] . In population biology, competition 
often has the range shorter than that of dispersal. In such a case, the mentioned 
homogenization occurs at nonzero mortality m. For instance, let a + = ckIr and 
a_ = /3I r for some positive a, /?, r, and i?, and also R > r. Here I r (x) = 1 if |x| < r 
and I r (x) = otherwise. Then ^2HJ) holds if 
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